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1. Introduction 

Calogero-Sutherland modelpQ El El SI El El El El EJ is a one-dimensional quantum 
system of particles with two-body interaction inversely proportional to the square of 
distance. From theoretical point of view, Calogero-Sutherland model has attracted 
extensive attention in relation to fractional exclusion statistics [10] [HI [12], Jack 
polynomials [T3[ E], Tomonaga-Luttinger liquid[15], collective field theory [T6l [TTj and 
matrix models[18j. An intriguing property of Calogero-Sutherland model lies in the 
fact that exact expressions of two-point dynamical correlation functions have been 
obtained for whole range of time and space [HI [19l EQl EH [22l [23l EH ESI [27]. 
The expressions for dynamical correlation functions are much simpler than those of 
Bethe solvable models [28]. Through exact explicit expressions of spectral functions of 
Calogero-Sutherland model, dynamical properties of one-dimensional quantum systems 
have been discussed[29, 30J. 

In this paper, we consider spin generalization of Calogero-Sutherland model[3Tll3~2"[ 
[33], which is the model of particles with spin 1/2 as internal degrees of freedom. In 
the following we call this model "spin Calogero-Sutherland model" while the original 
model "scalar Calogero-Sutherland model" . Spin Calogero-Sutherland models are even 
more important than scalar model in the sense that (i) the spin Calogero-Sutherland 
models have Yangian symmetry [341 135] as an internal symmetry, (ii) it reduces to the 
Haldane-Shastry model[36, 37J in the limit of infinite coupling parameter [38| [39] and (iii) 
it realizes the Tomonaga-Luttinger liquid (of particles with spin internal symmetry) in 
the systems with finite number of particles in the simplest manner [32| HQ]. 

Exact results of dynamical correlation functions such as hole propagator [HI |4*2] . 
density correlation [4*31 14*4"] and spin correlation function [4"3l l4"4l |4"5| |4"6] have been 
obtained in spin Calogero-Sutherland model. However particle propagator of this model 
has not been derived. With knowledge of hole and particle propagators, we can arrive 
at full understanding of the dynamics and elementary excitations. From the hole and 
particle propagators, furthermore, we can construct 1-particle causal Green's function, 
which provides a nontrivial starting point of many-body perturbation theory. The 
derivation of exact results applicable both to hole and particle propagators, therefore, 
is highly required from the above point of view. We address this issue in the present 
paper. 

There are two types of the expression of the wave functions of the spin Calogero- 
Sutherland model. The first one is given by Jack polynomial with prescribed 
symmetry [471 148]. This polynomial is constructed by partially symmetrization or anti- 
symmetrization of non-symmetric Jack polynomials [4*91 150] . which are the simultaneous 
eigenfunctions of the integrals of motion of this model, that is Cherednik-Dunkl 
operators [5T1 152]. Jack polynomials with prescribed symmetry form an orthogonal basis 
of the Hilbert space with a specific spin configuration. The second one is related to the 
Yangian symmetry [341 [35J of spin Calogero-Sutherland model. An orthogonal basis of 
the Fock space including spin degrees of freedom with fixed particle number is called 
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Yangian Gelfand-Zetlin basis[53J EI] 

Dynamical correlation functions of spin Calogero-Sutherland model were calculated 
in two methods according to the two types of eigenfunctions. One way uses Jack 
polynomials with prescribed symmetry as eigenfunctions, and the dynamical correlation 
functions are calculated by using the relations derived from the non-symmetric Jack 
polynomials. With this method, hole propagator [HI H2] has been derived by one of 
the authors and his collaborator. Recently [55], it was found that dynamical density 
correlation of spin Calogero-Sutherland model can be derived with Jack polynomials 
with prescribed symmetry and the method used in [56] . However, no formulas necessary 
to particle propagator have been obtained in the theory of Jack polynomials with 
prescribed symmetry. 

The other way uses Yangian Gelfand-Zetlin basis as eigenfunctions together with 
the mapping from Yangian Gelfand-Zetlin basis to symmetric polynomials |43] . The 
relations of this polynomial necessary to calculate dynamical correlation function are 
derived from those of Macdonald polynomials |14j. Density correlation function and 
spin correlation function have been derived with this method [4~3| |4~4"| . We refer to the 
latter method as Uglov's method. For a decade, it has been unresolved issue whether 
the Uglov's method is applicable to the calculation of 1-particle Green's function. In 
Uglov's method, wavefunctions of mult i- component particles are mapped to those of 
single-component particles. It is not obvious how the field operator in 1-particle Green's 
function is transformed under this mapping, in contrast to the density operator or spin 
operator considered in [43] . 

The main purpose of this paper is to show that Uglov's method is applicable to 
derive both hole and particle propagators. More explicitly, we give the transformation 
of field operator under the Uglov's mapping from multi-component system to single- 
component one. As an application, hole propagator is calculated by using this 
method. With introducing renormalized momenta and spin variables of quasi-holes, the 
expression for hole propagator in finite-sized system becomes much simpler than that 
derived in [12]. We confirm that the expression in the thermodynamic limit recovers 
the earlier result [4~2] . Furthermore we discuss spectral function of hole propagator 
for a specific coupling constant (A = 1 in the notations we will introduce in the 
following chapters). At this coupling parameter, the hole propagator of spin Calogero- 
Sutherland model is equivalent to the dynamical colour correlation function [431 |4"6] of 
SU(3) Haldane-Shastry model[3T l 132 1 [57] as shown by Arikawa[58j. In the same way as 
hole propagator, particle propagator can be mapped to that of single-component model. 
However, it is more involved to take the thermodynamic limit of particle propagator 
than hole propagator. Therefore, we report the calculation of particle propagator and 
discussion of the corresponding spectral weight in a separate paper. 

To outline this paper, the basic properties of Spin Calogero-Sutherland model are 
shown in section [2j In section [31 the mapping of the field annihilation operator is 
considered. Using the result of section [3j hole propagator of spin Calogero-Sutherland 
model is derived in section [H The expression for hole propagator is rewritten in terms 
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of quasi-hole rapidities in section The properties of the spectral function for A = 1 
are discussed in section |6j 

2. Basic Properties 

In this section, we review the basic properties of spin Calogero- Sutherland model. 
2.1. Hamiltonian and eigenf unctions 

Spin Calogero-Sutherland model is a one-dimensional quantum model which consists of 
N particles with spin degrees of freedom, moving along the circle of perimeter L. Each 
pair of particles has an interaction of inverse square type potential. Hamiltonian is given 
by 

" d* 2? r 2 v A(A + P 8J ) 

tT^ + ^ sin 2 [tt^ U 

where Xj is the coordinate of i-th particle, and Py is the spin exchange operator for 
particles % and j. The Hamiltonian (CQ) has one parameter, an interaction parameter A, 
which controls whole physical properties of the system. Physically, A < is unrealistic 
due to the collapse of the particles. In the following of this paper, we consider particles 
with spin 1/2, and A is restricted to be non-negative integer for simplicity. 

Let the particles be bosons for odd A and fermions for even A, following to the 
earlier works on the hole propagator (JH H2], and the boundary condition is chosen to 
be periodic. In the ground state, N/2 is taken to be odd (even) when A is even (odd) to 
avoid the degeneracy of the ground state. Introducing new variables z = {z\, ■ ■ ■ , z^} 
with Zi = exp[27iiXi/L], the wave function ^/({xj}, {cij}) for spin 1/2 model is written 
as the product $(z, a)^o,N( z ) of Jastrow type wave function 

*o^)=n^" i)/a n^-%) A ( 2 ) 

i i<j 

and a function <&(z,o~) of complex spatial coordinates z and spin variables a = 
{o"i, • • • , o~n}- The spin coordinate <jj takes 1 (2) for spin up (spin down). 

When A is even, particles are fermions and the Jastrow wave function ^o,n(z) is 
symmetric with respect to interchange of z^ and Zj. $>(z, a) then obeys the fermionic 
Fock condition, 

$(••■, Zi, a h ■ • ■ , zj, o-j, ■■■) = -$(■ • • , Zj, o-j, ■ ■ ■ , Zi, Oi, ■ ■ ■)• (3) 

When A is odd, on the other hand, the particles are bosons and the Jastrow wave function 
^o,n{ z ) is anti-symmetric with respect to interchange of Zi and Zj. Thus Q(z,a) obeys 
the fermionic Fock condition ([3]). 

The Jastrow wave function ^o,at is periodic under the translation x j — > Xi + L when 
A is even or iV is odd and anti-periodic otherwise. It thus follows that 

_ , . . periodic when A is even or N is odd , . s 

<&{z,a) is < . , (4) 

anti-periodic otherwise. 
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A basis for wave function satisfying Q and (jl]) is given by the Slater determinant of 
free fermions with spin 1/2 

' N 



u K:a = Asym 



i=i 



(5) 



for a set of momenta ft = (ki, • • • , kn) and spin configuration a = («i, • • • , ajv). Here 
ai — 1,2 means z component of spin of i-th particle being +1/2 and —1/2, respectively. 
Here one-particle spin function ip ai (o~i) is given by #3/2-a;,ov The symbol Asym[- • •] 
means anti-symmetrization of the function of z and a 

Asym f(z 1 ,a 1 ,---,z N ,a N ) 

= ^ ( _1 ) P f( Z P(l)' °>(1)> ' ' ' ' Z P(N), CP(N)), (6) 

Pes N 

where (— 1) P denotes the sign of the permutation P in the symmetric group Sn- 

When the basis function (j3j) obeys the periodic boundary condition, the set of 
momenta ft belongs to 

£n,2 = {k = (ki, « 2 , • • • , k n ) e £ N \ Vs g Z, | = s} < 2} , (7) 

which is a subset of 

C N = [k = («i, • - • , Kjy) G Z*|k< > ft m for i G [l,iV- 1]} . (8) 

We call the elements of by shifted partitions. When the basis function flSJ) obeys the 
anti-periodic boundary condition, ft belongs to C N2 , which is defined by 

C' N>2 = {k I k + 1/2 = (« x + 1/2, • • • , k n + 1/2) G £^,2} • (9) 

For spin 1/2 system, the definition ([7]) comes from the fact that a single orbital 
state can accommodate at most two particles. Furthermore a pair of particles with the 
same momentum ft, = ft i+1 cannot have the same spin state, i.e., a, 7^ c^ + i owing to 
Pauli exclusion principle. For a given set of momenta ft G £jv,2 or C' N2 , therefore, each 
spin configuration is specified by the element of W K defined as 

W K = {a = («i, • • • , a N ) G [1, 2]^! a { < a i+1 if ^ = . (10) 

For N = 2, W K is given by {(1, 1), (1, 2), (2, 1), (2, 2)} when k x > ft 2 , and {(1, 2)} when 
fti = ft 2 . Thus the basis function u KtCe is uniquely specified by (ft, a) G (Cn,2, W k ) 
under the periodic boundary condition and (ft, a) G (£jy- 2 , W*) under the anti-periodic 
boundary condition. 

Now we define the ordering between the basis functions. First we introduce 
dominance partial order [Tj 



r r 

v > n \v\ = \n\ and Vr > , > //j (11) 

i=l i=l 

between u, n G £jv,2 or e £jv2- Next we define the order for spin configurations as 

v v 
a; > a' <^ = a^, 

i=i i=i 

and nonzero ct i — a, at the least i is positive. (12) 
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For example, spin configurations with N = 3 and a i = 4 are arranged as 

(2,1,1) < (1,2,1) < (1,1,2). (13) 
The order of (k, a) is then defined by 

(k, a) > (k', a') k > k' , or k = k' and a > a'. (14) 

Uglov showed [43] that the excited parts of the eigenfunctions for the Hamiltonian (Op) 
is characterized by (k, a) as Q Kta (z,a), which can be uniquely defined by the following 
two conditions: 

(i) § KtCe (z, a) is expanded by u K ^a' satisfying (n',a') < (k, a) 

(/t',a')(<(«,a)) 

(ii) Orthogonal with respect to the norm (• ■ -)n,x 

^^)n,x = for {k\ a) ^ (k, a), (16) 
where the scalar product (• • -)n,x is defined by a weighted integral 



J_ fN 



(17) 



($(2,0") means complex conjugate of $(z, er)). The scalar product (jTTI) comes from the 
usual one: 

= S 2 f dx!-- /" i dx^'({x J },{a l })^({x,},{ ( x l }). (18) 

(71 =±1/2 (7JV=±l/2^ ^ 

The relation 

= ML JV ($',$) JViA (19) 

holds when \f = $^o,at and ^' = ^'^o.Af- The eigenenergy of the eigenfunction $ Kt0 ,^o,N 
is given by 

N 

E n (k) = (2tt/L) 2 ^ («< + A(iV + 1 - 2z)/2) 2 . (20) 

i 

For example, the ground state of iV-particle system is specified by [54] k°, a° with 

V4 2' 4 2' 4 2' ' 4 + 2' 4 + 2 J [ ] 

and 

a° = (l,2,l,---,l,2). 
The ground state energy E N (g) is then given by [2D E3] 

Ejv(g) = (2tt/L) 2 [(1 + 2A) 2 (M 2 - l)M/3 + A 2 M/2] (22) 
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with M = N/2. 

The above properties f|T5|) . f|T6l) and ffT7|) are used for the mapping of the 
eigenfunctions of spin 1/2 Calogero-Sutherland model in Uglov's method|43j. The 
polynomials $ KjQ (z, a) defined by the above two conditions (fT5|) and f|T6|) are known 
to coincide with Yangian Gelfand-Zetlin basis |431 131] . while Yangian Gelfand-Zetlin 
basis for spin 1/2 system is originally defined as the simultaneous eigenfunctions of 
the quantum determinants of Yangian Y(qIi) and Y(g[ 2 )[l3] [53], [31]. In the following, 
therefore, we refer to $ KjQ (z, a) as Yangian Gelfand-Zetlin basis. 



2.2. Macdonald polynomials and gh- Jack polynomials 

Macdonald polynomials [H] are symmetric polynomials with two parameters. As a limit 
of Macdonald polynomials, symmetric Jack polynomials and Schur polynomials can 
be derived. A lot of mathematical relations of Macdonald polynomial are known and 
provide important formulas of Jack polynomials, which have been utilized in calculation 
of correlation functions of scalar Calogero-Sutherland model[22l [23], [25] EH 127] . 

Macdonald polynomials themselves are eigenfunctions of Ruijsenaars-Schneider 
mo del [391 160] . which is a relativistic generalization of scalar Calogero-Sutherland 
model. The dynamical correlation functions of Ruijsenaars-Schneider model have been 
calculated with use of the properties of Macdonald polynomials[61J. 

As an index which specifies each symmetric polynomial, we define partitions as 
the set of non-negative integers arranged in non-increasing order. We denote the set of 
partitions with length equal to or shorter than N by 

K N = {v = {v u u 2 , ■ ■ ■ , v N ) G Z N \v x > v 2 > ■ ■ ■ > v N > }. (23) 

The monomial symmetric polynomial m v with a partition v £ A^ is defined by 
symmetrization of a monomial z v = z^z^ 2 ■ ■ ■ z v ^ as 

m v = Y^z{ w z^---zT\ (24) 

17 

where the sum is taken for all distinct permutations of the elements of v. 

Macdonald polynomial P u (z;q,t) for v £ is uniquely defined by the following 
two conditions [14]: 

(i) P u (z; q,t) is expanded by m M satisfying fi < v 

P v {z\ q,t)=m v + ^ ^"V ( 25 ) 

ju(<v) 

(ii) Orthogonal with respect to the norm (• ■ -)N,q,t 

(P^(z;q,t),P v (z;q,t)) Ntqt t = for \i ^ u, (26) 
where the norm of f[26l) is defined by a weighted integral through the function (x; g)^ = 

nr=o(i-^ r ) 



(/) 9)N,q,t — -j-jj 



' N 

n 

i=l 



dzj 

2-RlZi 



Y[Mf2l^f { z)g(z). (27) 
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Since symmetric Jack polynomial is the limit t = q x , q — > 1 of Macdonald polynomial, 
the above two conditions also define symmetric Jack polynomial uniquely by taking the 
limit t = q x , q — > 1 of the norm ( 1271) . Schur symmetric polynomials s u (z) for v G An 

8 ( z ) = Asym If* £? Z JL± (28) 

can also be obtained as the limit t = q — > 1 of P„(^; t). Uglov utilized the properties 
of Macdonald polynomials to calculate the dynamical correlation functions by mapping 
Yangian Gelfand-Zetlin basis to symmetric polynomials that are a limit of Macdonald 
polynomials [43J. These new polynomials are called gl 2 -Jack polynomials |43| 144] . and 
defined as 

Pj 2A+1 )(z) = lim P v {z). (29) 

q=— p,t=— p AA+L ,p— +1 

In the following, we call the limit in ([29]) "Uglov limit". From (1251). (|26l) and fl29|) . it 
follows that 

(i) Pj) 2X+1 \z) is expanded by m M satisfying /j, < u 

pV+ 1 )(z)=m v + <W (30) 

fi(<u) 

(ii) Orthogonal with respect to the scalar product {• • -}n,x 

{ p(2A + l) )P ( 2 A + l) W = ten^V. (31) 

The scalar product in (131]) is defined as 

{f^} N ,x = 
1 

M 



n 



2nizi 



A+l i \ A 



.1=1 

which comes from a limit of (|27|) . 

The two properties (!30|) and (131]) can be regarded as the defining properties of gl 2 -Jack 
polynomials. Alternatively, we can define gl 2 -Jack polynomials by 

P^\z) = s u + Cu,s, (33) 

H{<v) 

and (13T|) because (I33I) and (1301) are equivalent as shown below. 

Schur polynomials, which is a limit of Macdonald polynomials, can be expanded by 
monomial symmetric polynomials, and conversely monomial symmetric polynomials are 
written in the form of 

m u {z) =s u + 22 a uu'S u ', (34) 

v'(<v) 

from which 

P^\z)=m u+ J2 
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:m,o,-2) 



a = 

(1,2,2,2) 



a — 2k = 
(-1,0,2,6) 



v 



(3,0,-1,-1) 



Figure 1. An example of the transformation (|37| for N = 4. While (k, a) £ 
(£jv,2, W K ) or (C' N 2 , W K ), a — 2k is a shifted partition of the spinless fermionic system 
in reverse order, and v is a shifted partition of the spinless bosonic system. 



!/'(<!/) ft(<v) \ 

= + ^ C^s^ (35) 

f*(<I/) 

follows. 

So far Macdonald polynomials P„ and gi2-Jack polynomials pj 2A+1 - ) have been 
defined for a partition v G A^. However, it is convenient to extend the definition 
of P v and pj 2A+1 ) for v that belongs to Cn (or £jy). When i/ G £at is given by 

v = // - J = {fa - J, • • • , /i N - J) 

with an integer (half integer) J and a partition fi <E A N , we define P u and pj 2A+1 ) as 

P v (z) = (Z! • • •z JV )- 7 P /1 (z), Pj 2A+1) W = (*i • • •^)- 7 i'i 2A+1) W, (36) 

respectively. 



Mapping from Yangian Gelfand-Zetlin basis to glq,- Jack polynomials 

Uglov defined [33] the linear mapping Q between the set of functions spanned by with 
(k, a) G (£>n,2,W k ) or (C' N2 ,W K ) and the set of symmetric functions as fl(u Kja ) = s u , 
where the relation between (k, a) and v is given by 

Pi = a N+1 _i - 2k n+1 „ { -N + i (37) 

for the system of N particles with spin 1/2. An example of this transformation is drawn 
in Figure [U The properties of Q are listed as follows: 

(i) Isometry. The scalar product is preserved under the mapping Q. For functions 
&(z, a) and §(z,a), the relation 

(<z>',<s>) N>x = {n(<s>'),n(®)} N ,x (38) 

holds. 

(ii) For any symmetric function f(z) [33], the relation 

• • • , z N ) u K}a (z, a)) = f(z^ 2 , z^ 2 ) Q(u K>a (z, a)) (39) 

holds. 
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(iii) The correspondence between Yangian Gelfand-Zetlin basis and gl2- Jack polynomials 

n($«, a ) = pj 2A+1) . (40) 

Taking the mapping Q for both sides of (1151) and ( |T6l) in the conditions that specify 
Yangian Gelfand-Zetlin basis, there appear the defining relations of g^-Jack polynomial 
Py 2X+l \z) (|33|) and (I3TI) . and consequently, the property (iii) follows. The mapping 
Q can be interpreted as a transformation from a multi-component system to a single- 
component system. 

3. Transformation of the Field Annihilation Operator by the Mapping 

In this section, we consider the mapping of the field annihilation operator. In calculating 
density correlation function and spin correlation function, spin operator and density 
operator are expressed as power sum polynomials, that is, c-numbers|43j. Annihilation 
and creation operators in 1-particle Green's function, however, are not the case. We 
first make sure that the sum of the annihilation operators of spin f and j on the multi- 
component model are mapped to the annihilation operator on single-component model. 

3.1. Action of annihilation operator 

Generally, the action of the field annihilation operator on a wave function is implemented 
by fixing the coordinate of one of the particles in the wave function to that of the 
annihilation operator as 

ijj{x)^{x u ■ • • , xjv_i, x N ) = y/NC^Vfai, ■ • • , Zjv-i, x) (41) 

for single- component model, where ( = 1 for boson and ( = — 1 for fermion. Moreover 
for multi-component model, the action of the field annihilation operator with spatial 
coordinate x and spin coordinate o~(= 1/2,-1/2) is given by 

^ CT (x)^(xi, 0i, • • ■ , ZjV-l, 0JV-1, %N, Cr N ) 

= VN( N - 1 ^(x 1 ,a ir --,x N - 1 ,a N . 1 ,x,a). (42) 

In the following, we write ip-\(x) (ipi(x)) instead of ^1/2 (#) ("0-1/2 (^)) f° r notational 
convenience. First, we consider the action of (0|(O, 0) +0^(0,0)) on the wave function 
of spin Calogero-Sutherland model 

*(x,<T) = $(z,<T)# 0tN (z), (43) 

where \E , ,A r (-2) is the Jastrow type wave function defined in (j2J). We can take the action 
of one of the annihilation operators by restricting the spin configuration of the Yangian 
Gelfand-Zetlin basis of both sides of annihilation operator when calculating the matrix 
element. The sum of the field operators acts on the wave function ( 1431) as 

^ T (o,o)+^(o,o))* 

= C iY_1 v / iV : <&(zi,eri, erjv_i,Zjv = l,a N ) 

a N = l,2 
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x ^o,n( z i, * * ' > z N-li z n = 1) 
= ( N ~ 1 y/N ^2 ^(z 1 ,a 1 ,- ■ ■ ,z N ^,a N ^,z N = l,a N ) 

a N = l,2 

JV-1 

X 

1=1 

Next, we consider the action of the field annihilation operator on the wave function 



J I z't X/2 (zi - l) A ^ ,7v-i(^i, ■ ■ • , z N - X ). (44) 



f{z)*oA*) (45) 
where f(z) is a symmetric function of z and \l/o./v is given by 



N 



*w(z) = n^ (2A+lKJV - 1)/2 - ^) A+1 (^ + (46) 

i=l i<j 

This function is the Uglov limit (q = —p,t = —p 2X+1 ,p — > 1) of the ground state wave 
function for Ruijsenaars-Schneider mo del [591 160] . Acting the field annihilation operator 
on (j4"5]) . we obtain 

^(0,0)/(z)* o ,w(s) 

= C^ViV/C^i, ■ • ^Zn-uZn = l)^ ,7v(^i, • ■ -,z N -i,z N = 1) 
= £ N - 1 VNf(z 1 ,-- -,z N - 1 ,z N = l) 



N-l 
i=l 



II ^ A_1/2 (^ - !) A+1 (^ + 1) A *o,jv-i(«i, ■ • • , Ziv-i). (47) 



We introduce ip s (0, 0) and -0(0, 0) through the relations: 

&(0, 0)$ = (* 0> iv-i) _1 V.(0, 0) $^ ,v, s =T or | (48) 
$(0, 0)/ = (^iv-i)" 1 ^(0, 0) f4? , N . (49) 
In terms of (1481) and (1491) . fj44p and (1471) are rewritten as 
(^ T (0,0) +^(0,0))$ 

JV-1 

= C iY ~ 1 v / iV ^ $( 2l , r 1 ,---, 2iV _ 1 ,a w _ 1 ,^ = l ! a Ar ) Hz; X/2 ( Zl -l) x 

o-jv=1,2 i=l 

(50) 



and 



^(0,0)/(z) 



JV-l 



C^VNfizu ■ • • , Ziv-i, ^ = 1) I] ^ A " 1/2 (2< - l) A+1 (^ + 1)\ (51) 



1=1 



respectively. 
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3.2. Transformation of the field annihilation operator 
In this subsection, we show that 

n ($ T (o, o) + ^(o, o))$) = (_i)<"-i>* (]J z \ /2 ] ^(o, o)n(*). (52) 

To derive fl52l) . it suffices to show 



i=i 



'N-l 



n (^(o,o) + ^(o,o)K J , 



(53) 



i=l 



with ( 1371) . This is because (i) Q is a linear operator, (ii) u Kj0l is a basis of $ and (iii) 

Since m KjQ , is a Slater determinant of N free fermions with spin 1/2, u K>a with one 
of the coordinate fixed can be expanded by Slater determinants of N — 1 free fermions 
with spin 1/2 as 

22 u *A z ii °i 3 • ■ • > z n = i, o") 



(7=1,2 



E 

<T=1,2 



Vai(0"l) 
^Ar-lVaiiCjV-l, 



A' 



^(-l) l M... jKi _ liKi+li .. v .. iQ8 _ liQi+li .../s(zi, (Ji, • • • , Zjv-l, CTjv_i). 



(54) 



i=i 



Therefore, 



n((V T (o,o) + Vi(o,o)K 



AT-1 



% A/2 U; - d a 



ctjv=1,2 



1=1 



AT-l 



A? 



C^Vn I) 



i) A £(-i) 



i=l 



i=l 



^ ,K i _i,«i + i ) - > -,a i _i,ct i+1 ,-(^l 5 CTl, • ■ ■ > %N-1, O'N-l)) 



N-l 



N 



c 



N-l 



1=1 



i=i 



X S i/l+l,---,l/]V-i + l^JV_i + 2,---^]v( 2 l' * " J ^JV-l)- (55) 

Here we use the property flHUl) of f2. We note that though v is defined by v± = 
ajv+i-j — 2/t7v +1 _j — N + i ioi iV-particle system, the partition for Schur polynomial 
on the right-hand side of (|55|) has N — 1 elements. 

Since Schur polynomial is defined by Slater determinant of spinless fermion, the 
right-hand side of fl55|) can be described by a Schur polynomial of N variables with one 



Derivation of Green's Function of Spin Calogero- Sutherland Model 
of the variables fixed as 

s u ( z i, ■ ■ ■ , zjsr-i, z N = 1) 



13 



Asym 



vi+N-l i^jv-i+l-i v N 

Z\ ■ ■ ■ z N _ 1 1 



Asym \jz\ 



w-i . „1 10 



2 



ui+N-1 
1 



2f 



iV-1 



-N-l 



N-l 

in*- 1 ) n 1 

l<i<i<iV— 1 



-I -1 



Z% Zj t 



1=1 



Eii(-l) J Asym 



■ z, 



Vi^+N-ii-l) v i+1 +N-(i+l) _ 1 Ar _ 1 



Asym 



N-2 



J 



"JV-1J i=i 
N N-l 

= / J(~l) ts ^i+l,-^jy-»+l,^-i+2,-,^w( z l; ' " " j z iV-l) JJ (z» - 1) 1 - (56) 
i=l 

Thus the right-hand side of fl55|) is rewritten as 
r.h.s. of (55) 



i=l 



N-l 



i=l 



'N-l 



(57) 



i=l 



From (155]) and (15T|) . the relation (1531) follows. Using (152]) . the matrix element of field 
annihilation operator is obtained as 

$ T (0, 0) + ^(0, O))^)^^ 

= {fi^aO, ^((^t(0, 0) + ^(0, 0))® K , a )} N -l,X 



'N-l 



_ 1){ N-1)X p(2A + l) ; J-J z l/2 ^ (0)0 )P( 2 ^ 1 ) 



i=l 



AT-1,A 



(-i)^-^{p^^(o,o)Pi 2 ^ 



AT-l.A 



(5f 



Here we note that the relation between v' and («', a') is defined by ( }3~Tj) with the number 
of particles being A^ — 1 while the relation between v and («, a) is that with A". 



4. Combinatorial description of hole propagator 

Hole propagator is one of the 1-particle Green's function defined as 
t , (g,iV|^(rc,i)*l(0,0)|g,iV) 

Gfet) = feM ' (M) 

Here |g, N) is the state vector of A^-particle ground state, whose wave function is given 
by a°^o,N- The scalar product in (1511 is the conventional one (see (fT81) ). 
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We rewrite (1591) in terms of gVJack polynomial. First of all, a complete set of the 
state vectors is inserted between the creation operator and the annihilation operator in 
the numerator of the right-hand side of ( 1591) . Next the complex conjugate is taken so as 
to alter the creation operator to the annihilation operator. Third, we use the relation 

((K,a),JV-l|^(0,0)|g,iV) = ((«:,«), iV-l|(^ T (0,0)+^(0,0))|g,iV)(60) 

when z component of the total spin S z of the state («, a) is larger than that of the ground 
state by 1/2. The symbol \(k, a),N — 1) denotes the state vector of (N — l)-particle 
state whose wavefunction is $ KjQ \l/o,v-i- The scalar product between the state vectors 
is then represented by the scalar product between functions of Yangian Gelfand-Zetlin 
basis. Finally we act the mapping Q on the wave function and obtain 

G~(x,t) 

-^t +i P KX I ((*, <*),N-1 1 (^(0, 0) + ^ (0, 0)) |g, iV) | 2 



e 



((K,a),N-l\( K ,a),N-l)-(g,N\g,N) 

Stot=+V2 



i e _ i ^ +iPK j($ K , a ,(^(o,o) + ^ T (o,o))$g)^ 1 , A 



2 



$g , $g) jv,a 
=-*-Stot=+ 1 /2 



LN ^ rp(2A+l) p (2A+l) x rp(2A+l) p (2A+l)-, ' 

s ' t -S t Z ot =+l/2 

(61) 

The variables of summation v is related with (k, a) via 

Ui = a N -i-2K N -i-N + l-i, ze[l,JV-l]. (62) 

The energy ui K and momentum P K are described in terms of k G Cn,2 or C N 2 as 

/ u R = E n ^(k) - E N (g) 

2vr v (63) 

K ~ T / J 1 



L 

i 



with 



B -W = (T)'?(- + M f M )' (64 » 

and (J22D- n maps $ K<V o to F$ A+1) = P^ X+1) with 

n f N N N \ 

" =(" + 2 --T + 2 -"'--T + 2 J- (65) 

and we obtain gl 2 -Jack polynomial of the ground state as Pg X+l \z) = Y\ i z^ N ^ 2+2 . The 
restriction on the sum is considered later. 
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k = (5,5,4,4,4,3,1,1) <£> 
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Figure 2. Partition and Young diagram. 



^.1. Expansion by gh-Jack polynomials 

The matrix elements of correlation functions in the Sutherland model are expressed in 
terms of partitions. Partitions can be expressed graphically by Young diagrams, whose 
correspondence with partition is drawn in Figure [3 From the top, k% squares are placed 
in the first row, and k 2 squares in the second row, and so on. Each square is specified 
by two-dimensional coordinates, labeling the square at the upper left by s — (1, 1). The 
first coordinate indicates the vertical position and the second the horizontal position. 
The length of a partition is defined by the number of nonzero elements in the partition, 
and equal to the length of the first column of the Young diagram. 

Some relations of symmetric polynomials used in calculating correlation 
functions are described in terms of the variables defined on the Young diagram, 
a(s),a'(s),l(s),l'(s) for s = as 

a(s) = Ki - j, l(s) = Kj - i, 
a'{s)=j-l, l'{s)=i-l, 

where Ki is the i-th element of the partition k, and «;'• is the length of j-th column in 
the Young diagram (Figure [2]). 

In calculating the numerator in the right-hand side of ( l6Ti) . the term — 
l) x+1 (zi + 1) A appears in ^(0, 0)P^ X+1 \z) and it can be expanded by gl 2 -Jack 
polynomials. The expansion formula is obtained from the corresponding relation of 
Macdonald polynomials [14j 

n^= E «^: 5 , t ), (66, 



where 



i (rt t] = n ( tns) - ^ 

h»>(t,q)= H (l-g-W+^W) 



(67) 
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Figure 3. C 2 {^) and H 2 {^)- The white squares belong to G^ifj) and H 2 {p) 
respectively. 



Taking the limit q = —p , t = —p IX+1 , p — > 1 for f[66|) . we obtain 

N-l 

n (i - + Zi t = Yi p r +i) & ■ (-i) iH+Ei,(s) 



i=l 



m£A jv _ 1 
s.t. |o 2 (m)I+|h 2 MI=|mI 



J] (a'(s)-(2X + lW(s) + l)) 



(68) 



Yl (a(s) + 1 + (2A + l)Z(s)) 

where C 2 (fi) , H 2 (fi) are the subsets of the squares in Young diagram D(fi) defined 
respectively as (Figure [3]) 

C 2 (^) = { s E D(fi)\a\s) + l'(s) = mod2} (69) 
H 2 (ii) = {s e D(p;)\a(s)+l(s) + 1 = mod2}. (70) 

The set A \ B means the complementary set of B in A. In the right-hand side of 
(16"8"|) . the sum of /x is restricted to the partition satisfying |C 2 (/i)| + \H 2 {ji)\ = The 
numerator of (f6Tj) is written in terms of the variables of Young diagram and the norm 
of intermediate states as 



» £ 



<V-l/2,^-(iV/2+A-2)-l/2 

S.t. |C 2 (m)| + |H 2 (m)| = W 



(71) 



x 



f n (a'(s)-(2X+l)(l'(s) + l))\ 

seD(ti)/c 2 (iJ,) 

Y[ (a(s) + l + (2A + l)Z(s)) 



/ 



r p (2A+l) p(2A+m 

I A* ' M JiV-1, 



(72) 



The norm of gl 2 -Jack polynomial is also obtained by that of Macdonald polynomial [1 

jp(2A+l)p(2A+l)j 



TV, A 



„(2A+1,2) 
-N 



n 



a'{s) + {2\ + l){N -l'{s)) 
a'(s) + l + (2A+l)(jV-Z'(s)-l) 
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n a(a) + 1 + (2A+ l)f(a) 

JJ rt aW + (2A + 1)( ' w + 1) ' 

where c N is the norm of the ground state of the system of N particles with spin 
degrees of freedom 2 

,(2A+1,2)_ "Q C (2X+l) {j _ i)j (74) 
l<i<j<N 

(2A+l)(fc + l) W (2A+l)(fc-l) ,\ 

■ 2 - ^ \ fc = lmod2 

r((24ili+i) 2 

(75) 

k = mod 2. 



C 7V 



C( 2A+1 )(A;) = < 



p/ (2A+l)(fc+l) 1 W (2A+l)(fc-l) 1' 



r (lH*+i2*)r(lH*+i2*+i) 

To summarize this section, hole propagator is written in terms of the variables defined 
on Young diagram as 

_ r((A + i/2)iv-A)r(A + i) 

1 ' j T((A + l/2)iV)L 

X V p-^t+iP^x X l Y ^) (7(]) 

with a = 1/(2(2A + 1)). Cj^ and P M are, respectively, given by tD M = c<j k and P M = P K in 
f l63|) through the relation /i, = a^v-i — 2kn^ — iV/2 + A — 1 — i. The prefactor in front 
of the summation comes from c^^ 1 / (c^ X+ L) . We have introduced the following 
notations: 

X M = J] (- aa '(s) + (l'(s) + l)/2), (77) 

s eD( M )\c 2 ( At ) 

F M (r)= J] ( aa '(s) + r + (N-l-l'(s))/2), (78) 

s6C 2 (m) 

^(r)= J] (aa(s)+r + l(s)/2). (79) 
4-2. Restrictions on \i 

Combining the condition on v with that on /i, that is, 

(i) z component of total spin of the intermediate state is larger than that of the ground 
state added by 1/2, 

(ii) |c 2 (//)i + i# 2 ( M )| = M, 

there are two conditions for \i to satisfy in the sum. On the other hand, the factor 
a'(s) — (2A + l)(Z'(s) + 1) in the product of (1761) implies that the partitions including 
the square at (1, 2A + 2) have no contribution in the sum of (IToT) . Hence in the following 
we consider only the partition with [L\ < 2A + 1 and length equal to or shorter than 
N — 1. An example is shown for A = 2 in Figure HI We define the subsets of 2 A + 1 
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Figure 4. Young diagram of a quasi-hole state for A = 
= {8,6,4,3,1}, P = {1,3,4} and Q = {2,5}. 
to C2(n) and shaded one belong to D(fi) \ 



2. In this diagram, 
White squares belong 



columns of the Young diagram by P and Q as follows. The jth column belongs to P 
when the difference between j and the length of the jth column ylj is odd, and Q is the 
complementary set of P 

P = {jG[l,2A + l]|^.-j:odd} (80) 

Q = {je[l,2A + l]|/4-j:even}. (81) 

In Figure H = {8,6,4,3,1} and hence P = {1,3,4} and Q = {2,5}. 

We define n as the number of the columns in the set Q. Then |C*2(/i)|, |i?2(^)| a re 
written in terms of n as 

\CM\ = \tt±±Jl (g2) 

= l^|-n-A(2A + l) + 2(2A.n- n ( n -l)) ^ 



These relations are derived in Appendix A From fl82l) and fl83l) . the condition 



|C 2 (/x) | + l-f^Z-OI — |a*| can be described as the condition on n by 

\C 2 (ji)\ + \H 2 (ji)\ = \(jl\ & n = X,X + l. (84) 
Next, z component of the total spin for a partition [i is written as|43j 

s» = -\cm\ + \dm\cm\ + ~ 

= \-n+ 1 -. (85) 
Therefore the condition for total spin is also written in terms of n by 

= 1/2 <£> n = A. (86) 
The relation fl85l) gives the meaning of the relation (|84p that intermediate states fi that 
arise by acting (ipi + V't) on ^he wave function of the ground state contain only the 
states with = +1/2 or —1/2. Then the relation (IHBT) imposes the restriction to take 
only one of the annihilation operator. Putting together all conditions, the sum of fi is 
taken over the partitions satisfying 

Hi < 2A + 1 and length equal to or shorter than — 1 (87) 
n = A. (88) 
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5. Quasi-hole Description of Hole Propagator 

In this section, we rewrite the expression for hole propagator in terms of rapidities and 
spins of quasi-holes. 

5.1. Rapidities and spins of quasi-holes 

From ( 1H71) . N — 1 particle states relevant to the hole propagator ( 1761) are parameterized 
by the set of length of each column ■ ■ ■ , /4a+i} an d the set of "spin variables" 

{a 1 , ■ ■ ■ , CT2A+1} , 

the entry of which is defined by 

for i G [1, 2A + 1]. For later convenience, we introduce auxiliary notations n' = N — 1, 
/4a+2 = and 

(T = 1/2, <7 2A+2 = -1/2. (90) 

Regarding the length fi' of "0-th" column as iV — 1, which is odd and that of 2 A + 2-th 
as 0, the definition (1901) is a natural extension of ( 1891 ). 

Furthermore, we introduce the renormalized momentum 

* = "i — + ^xTT (91) 

for j G [0,2A + 2]. In terms of (f59"j) and (IDT!) for j G [1,2A + 1], excitation energy of 
eigenstates relevant to hole propagator is written in the form of free particles. Matrix 
element appearing in hole propagator is written in terms of {fij, aj} for j G [0, 2 A + 2]. 
We will see in the following that jl'- + aj and aj can be interpreted as the rapidity and 
spin of j-th quasi-hole, respectively. 

5.2. Hole propagator in finite-sized systems 

The excitation energy uj^ = uj k and the momentum P M = P K are described by 

My , ^ 2 ^7 ^ + ^^ 2 , 4vr 2 A(A + l) 
^ = -(2A + l)^l 3 - J + — (92) 

j=i \ / 

N-l 2A+1 

= (27T/L) 53 k, = - j3(7r(/i; + a 3 )lL). (93) 
i=i i=i 



The derivation of ( 1921) and ( 1931) is given in Appendix B The matrix element in the 
hole propagator (!76|) can be described in terms of the renormalized momenta and spin 
variables as the energy spectrum. 



Derivation of Green's Function of Spin Calogero- Sutherland Model 20 
X^, Y^ir) and Z^{r) defined in (1771) . (1T8"j) and (1?9~j) . respectively, are described as 

y rb7(2A + i)] ' m 

Y (r) = TT r[a + r + iV/2 + (j-l)/(2A + l)] 

W 1 = 1 r[(ft - # + 1 + <W,)/2 - a + r] ' lJ " j 
^(r) = (r[r + l/2])-( 2A+1 ) 

2A+1 

x II r K$ - - ^a 2X+2 )/2 + r - a + 1/2] 

r[(g -^-5 )/ 2 + r-a + l/2] 



A A IW- )/2 + rl 

l<j<fc<2A+l Lvr J j fe ' ' 



as shown in Appendix B With use of (1941) . ( 1951) and ( 1961) . the expression ( ITS"]) is rewritten 

as 

G~(x,t) = K x {N)d ^2 ^ ^(Ej ^).V 2 ex P[-^^ + 

0</4a+i<-</4<-N-1 

xF({/2;,^}), 



(97) 



with the particle density d = N/L and the constant K\(N) 

r((A + l/2)iV - A)r(A + 1) (r[((A + 1)/(2A + 1))]) 2A+1 



K X (N) 



J= 

and form factor 

2A+1 



r((A + i/2)iv - A) n?=t 1 (r[j/(2A + l)]) 2 

^ r[y7(2A + l) + iV/2-a] 

11 rbV(2A + 1) + N/2 - 1/2] 1 



j-j r[(&-/2; + Vo)/ 2 ] 



' rp' -^ + i + Vo)/2 



x 



r[(g - + 1 - 



n 



L r[(# - /}' 2A+2 + 1 - <V, CT2A+2 )/2 + 1/2 - a] 

rw. - # + 5 CT . a J/2 + a]rw. - ft + ^J/2 + 1/2] 



X 3 



n 

l<i<fc<2A 



1 r[(/2j - ft - * w )/2 + l/2]r[(ft - ft " <WJ/ 2 + 1 - a] ■ 

(99) 



Derivation of Green's Function of Spin Calogero- Sutherland Model 



21 



5.3. Thermodynamic limit 

Changing the variables to new ones by 

(100) 



II 5 1 - Uj 



N 2 ' 

which are finite in the thermodynamic limit. Using the relation T(N + a)/T(N + b) — > 
N a ~ h (N — > oo), ( 1761) converges to an expression with finite value. The result in the 
thermodynamic limit N — > oo{N/L = const.) is written as 



(2A+1 „i \ 

n e / °H 
4=1 ^il/2^ 1 / 

2A+1 2A+1 



>< n (i-^ 2 ) _a/(2a+i) n w-u^-^^) 

3=1 j<k 

/ , a2 2A+1 \ / 2A+1 

xexp ^(2A + 1)(^-J E u ?J* + < ^yE 



| X 



(101) 

with the overall constant 

K> T{X + 1} W r ((A + l)/(2A + l)) nn2l 

A 4(2A + l)*r(2A + 2) r(j/(2A + l)) 2 ' 1 J 

The exponents of interparticle part of the matrix element in (110 II) depend on whether 
spins of a pair of particles are parallel or anti-parallel. The result coincides with the 
earlier result 141, 421. 



6. Spectral Function for A = 1 

Spectral function for hole propagator is defined as a Fourier transformation of G~(x, t): 

-I /*oo /*oo 

A-(p,e) = — dx dte i ^ t - iPx G-{x } t). (103) 

J —oo J —oo 

We calculate the spectral function at A = 1, which is the easiest nontrivial case in this 
model. The spectral function has non-zero value at finite area in energy-momentum 
plane (Figure [5]). The area is enclosed by four parabolic lines, and each of them can 
be interpreted in terms of quasi-hole picture. The upper edge e = — p 2 + 9(7r<i) 2 /4 
with \p\ < 3nd/2 corresponds to the excited state with three quasi-holes having a same 
momentum. The three lower edges 

e = -3p 2 + 3(vrrf) 2 /4, \p\ < nd/2, 

e = _3(p + ird) 2 + 3(7rd) 2 /4, -3vrd/2 < p < -nd/2 

e = — 3(p - nd) 2 + 3(vrc/) 2 /4, nd/2 < p < 3vrrf/2 

correspond to the states that one quasi-hole is excited while the other two quasi-holes 
reside at Fermi points. 
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The intensity of the spectral function is also drawn in Figure |5j The spectral 
function diverges at two lower edges e = — 3(p±nd) 2 + 3(7id) 2 /4: for Tid/2 < \p\ < 3nd/2 
and two other parabolic lines e = — 3(p^fird/2) 2 /2 + 3(7rd) 2 /2, with —Tid/2 < p < 3ird/2 
for upper sign and —3itd/2 < p < —ird/2 for lower sign. At the lower edge 
e = — 3p 2 + 3(7T(i) 2 /4 for —ird/2 < p < 7cd/2, the spectral function becomes zero 
and arises as (e + 3p 2 - 3(7r<i) 2 /4) 1/3 , and at e = -3(p ± ird) 2 + 3(nd) 2 /A it diverges 
as (e + 3(p ± iid) 2 — 3(%d) 2 / '4) -1 / 3 . At the middle line, spectral function diverges as 
|e + 3(p =F nd/2) 2 /2 - 3(nd) 2 /2\~ 1/6 , which will be derived in |Appendix C| At upper 



edge, the spectral function takes finite value proportional to (1 — (2p/3ird) 2 )~ 1 . Arikawa 
pointed out [58] that the hole propagator of SU(2) Calogero-Sutherland model for A = 1 
is equivalent to the dynamical color correlation function of SU(3) Haldane-Shastry 
model[l5l HH1 [58]. Although the singularities of the upper edge and lower edges have 
been obtained in earlier works[46, 58J, our calculation shows that the spectral function 
at the middle line diverges with exponent 1/6, contrary to |46j . 



7. Conclusion 



To calculate the particle propagator of spin Calogero-Sutherland model by Uglov's 
method, we transformed the field annihilation operator on Yangian Gelfand-Zetlin basis 
by the mapping Q and proved that it becomes also the field annihilation operator on gl2- 
Jack polynomials. This ensures the possibility of calculating 1-particle Green's function 
using the Uglov's method. 

Next, with use of this method, we calculated hole propagator for non- negative 
integer interaction parameter by taking the restricted product over Young diagrams of 
intermediate states. Thermodynamic limit of hole propagator was also taken, and we 
confirmed that the result obtained here coincides with that of the former results with 
Jack polynomial with prescribed symmetry. Spectral function for A = 1 was calculated, 
and drawn in energy-momentum plane. There appear the divergences of intensity on 
one and two quasi-hole excitation lines. 
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The calculation of particle propagator will be published in a separate paper. 
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Appendix A. Restriction on the Product 

The number of the squares of the subset of the Young diagram, |C 2 (/x) | and |-ff 2 (//)|, 
can be uniquely described in terms of n, the number of the columns which belong to Q. 
This can be verified by describing |C 2 (/z) | — \D((i) \ C 2 (/x) | and \H 2 (fj,)\ — \D(/j,) \if 2 (/i)| 
in terms of n, where A \ B means the complement of B in A. 

In each column and each row of the Young diagram, the element of C^CaO and 
that of D(fi) \ C 2 (fi) are aligned alternately, and at the upper left of the diagram 
s = (1, 1) G C 2 (yu). From the definition of the set P and Q ( flHUl) and (jHTjl ). j-th column 
with j £ P has even number of squares for odd j and odd number for even j. Thus the 
square at the bottom of the j-th column in P is an element of D(fj,) \ C 2 (fi). For n = 0, 
all the columns belong to the set P, and therefore |C 2 (//)| — \D(li) \ C 2 (/i)| = —A. The 
change of a column in P to Q increases |C 2 (/i) — -D(/i)/C 2 (/z)| by one (Figure ED). We 
thus obtain 

\C 2 (fi)\ - \D(fi) \ C 2 (/i)| — —\ + n. (A.l) 

Next we describe \H 2 (n)\ — \D{jj) \ H 2 (fi)\ in terms of n. First we consider a 
partition having no adjacent columns with the same length, and then we consider the 
general case. In each column, the square at the bottom s = (fjfj,j) is an element of 
D(fi) \ if 2 (/i). The element of i? 2 (/i) and that of D(n) \ H 2 {ii) are aligned alternately 
except for several rows. The exceptions occur at the // fc -th row and fi' k + 1-th row in 
j-th column with k = [j + 1, j + 2, • • • , 2A + 1], where two elements of the same subset 
H 2 (n) or D(fi) \H 2 (fi) are aligned vertically (Figure lA"2l) . Thus we remove the fi'j + 1-th 
row with j = [2, 3, • • • , 2A + 1] from the original Young diagram, and make the diagrams 
D(ji) in which the element of H 2 (fi) and that of D{^i) \ H 2 (n) are aligned alternately 
without any exceptions in each column. The removed 2A rows make a new diagram 
D(Afi) with length 2A and (A//)i = 2A (Figure \M>- 

For n = 0, each column of D(jl) has even number of squares. Therefore, the 
contribution to \H 2 {n) \ — |-D(/i) \ H 2 (fi)\ from fx is given by 

|fr a (A*)UHa(/0\#2(/0U = o. (A.2) 

When a column in D(li) changes from P to Q, the number of columns in D(fi) with 
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n 







n > 



Figure Al. The relation between the number of the squares of Cz{p) and D{fi)\C^{ii) 
and n. White squares are the elements of C%{p), and shaded squares are that of 
_D(/i) \ C^Cm)- For n = 0, the bottom of each column are shaded squares (the elements 
of the -D(/i) \ C2(/i)), and for n > 0, there appears the white square (the elements of 
the (/•«)) in n columns. 



odd length increases by one. We thus obtain 

\H 2 {v)\»-\D{im)\HM\i = -n (A.3) 

for n > 0. 

Now we consider the contribution of D(A/j l ) to \H 2 (fi)\ — \D(fj,)\H 2 (fj J )\. The square 
s — [ji k) with 1 < j < 2A , 1 < k < 2A + 1— j in D(A/i) comes from s' = (/4a+2-j + 1j ^) 
in D(n). With use of a(s') = 2A + 2 — j — A; — 1 and Z(s') = — [i 2 \+2-j ~ 1> ^ follows 
that 

a(S)+/(S) + 1 -\ (*, 2A + 2 - j) e(P,Q) or (Q,P) ' (A ' 4) 

and hence s = (j, k) G #2(/-0 in D(Afi) <^> (A;, 2A + 2 - j) G (P, Q) or (Q, P). Since two 
coordinates k and 2A + 2 — j vary with 1 < k < 2A + 2 — j < 2 A + 1, 

|# 2 (»|A^ = n(2A + l-n) (A.5) 

|P>(//) \ F 2 (//)| Aai = n(n - l)/2 + (2A + 1 - n)(2A - n)/2 (A.6) 

and 

\H 2 {h)\ah ~ \D(fj) \ P 2 (/i)| AM = ~2n 2 + 2(2A + l)n - A(2A + 1) (A.7) 

holds. From ( 1A.3I) and ( 1A.7I) . we obtain 

|F 2 (//)| - \D(ji) \ H % (ji)\ = -2n 2 + (4A + l)n - A(2A + 1). (A.8) 

Next, we apply the result to a partition having two or more columns with the same 
length. When a partition have m columns with same length, we extract the 2[m/2] 
neighboring columns from the original partition, where [n] means the maximum integer 
that does not exceed n. Since the extracted 2[m/2] columns have the same number of the 
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1/2 



/4 



/4 



/i' 2 + l 



2A 



+ 



2A 



D(A/i) 



Figure A2. Young diagram D{fi) is decomposed to two diagrams D(A/i) and D(jl), 
White squares are the elements of #2 (/•*), and shaded squares are that of D(/i)\i?2(/u)- 
In each column of -D(A/i), the element of i?2(/-t) and that of -D(/i) \ -^(aO are aligned 
alternately. 



elements of H 2 (pi) and £>(/i) \_/7 2 (At), there is no contribution to |if 2 (/-0| — \ D(/i) \H 2 (p:)\. 
We can use the result (IA.8P for the left partition with 2A + 1 — 2[m/2] columns and obtain 
the equation by replacing A in (1A.8I) by A — [m/2]. The above discussion can be applied 
to a partition that has more than one set of columns with the same length. These 
applications, however, do not change the result on the condition |C 2 (/i)| + \H 2 (fi)\ = 

Appendix B. Energy Spectrum and Matrix Elements 



In this appendix, we derive (|92|) and (1931) in Appendix B.l and (|94|) . flMl) and (1961) in 



Appendix B.2 



Appendix B.l. Energy spectrum 

In the excitation energy u; K in fl63|) . the energy 

iV-l 



= (2tt/L) 2 ^ («i + \{N - 2z)/2f 



of AT — 1 particle state k is rewritten as 



JV-l 



= (tt/L) 2 (-^ + A + (1 + 2A) (i - AT/2) + aiV _, - l) 2 , (B.l) 



i=i 



with use of 



z/j = ajv-i - 2/tAr_j - N + l + i 



(B.2) 
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(which is (j37p with replacement of N by iV — 1) and 

l i = v + N/2 + \-2 (B.3) 
coming from ( 17T1) . In ( IB. II) . the value of ctN-i is given by 
= / 2, (i,fn) G C 2 (//) 

N ~ l \ 1, e DM/CM ' 

which follows from 

oiN-i = Hi — i + 2/€jv_j + N/2 + 1 — A 

even even 

= a\i, fii) + l'(i, /ij) mod 2. (B.5) 

For i G [yUi + 1, AT — 1], we regard (i, 0) as an element of C^/i) (D(fi) \ C^ip)) when z is 
even (odd). With use of (1B.4|) . E n _i(k) is rewritten as 

En-i(k) = ?r2(1 ^ 2A)2 (*W*) + ^_ x (k)) , (B.6) 

with 



i€[l,JV-l]; s.t. (i,/ij)eC 2 (/i) 



and 



i6[l,JV-l]; s.t. (i ift )eD( /1 )\C 2 (/ 1 ) V 7 

Now we consider It is convenient to decompose the sum with respect to i as 

2A+1 

E E < B - 9 > 

s .t. (,,, )e c 2W gt { J+ eC ^ (M) 

In the interval i G [1 + /i^], i's satisfying (i,fii) G C 2 (/i) align alternately. The 
minimum (maximum) value im in (j) (^ ax (i)) °f * m [l+A*j+i, /^] satisfying (z, //j) G C^/z) 
are given, respectively, by 

G*0") = + ^'+1 + 3 / 2 ' ^max(j) = H'i - ^ " 1/2, (B.10) 

as shown below. First we note that (i,fii) = in the interval i G [1 + 

When j + 1 G Q, — (j + 1) is even and hence (1 + fJ,j +1 ,j) G C 2 (/i) and 

When j + 1 G P, on the other hand, — (j + 1) is odd and 
hence (1 + j) G \ C 2 {n) and imin(i) = 2 + Thus we arrive at the first 

equation of flB~T0|) . We can obtain the second equation of (IB. 101) in a similar way. In 
term of i° in (j) and i° ax (j), £'n-i( k ) is written as 

N A — J + 1 

2~ + 2A + 1 



E E (i-T + 

3=0 i=P. fj) fl. (j)+2,- 

min w/ ' min w ' ' 
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2A+1 / x _• , -, \ 3 

6 



1 2A+1 / \ ' I 1 \ 3 

. 2A+1 

+ 6EH-x(j) + *°„(j)-2). (B.ll) 
In the second equality, we have used following formula: 

E (* + A ) 2 

i=m,m+2,m+4,--- ,n 

— - \(n + 1 + A) 3 — (m — 1 + A) 3 —n + m — 2] (B.12) 
6 

for n, m satisfying n — m being a positive and even integer. The expression ( IB. Ill) can 
be further rewritten as 

1 2A+1 

£ n-M = gE [(# - ^') 3 - + °* + V(2A + l)) 3 + 2a, 
i=i 

M(M 2 - 1) A 2 M 1 

+ 3 + (2ATiy + 6 (R13) 

with use of (IB. 10|) and (19T1) . We can rewrite the expression (IB. 81) for ^_ 1 (k) in a similar 
way. First ( IB. 81) is rewritten as 

2A+1 / , r . . < 



j=0 ie[l+Mj- +1 >^-], 

s.t. (i,n l )eD( M )\c 2 (M) 



2 2A + 1 



2A + 1 CaxO') / tit , ■ \ 2 

N X — j 



E E 

j=0 i=i m . (j),i m . (j)+2, 

mm \ J mm ^ * 1 ' 

2A+1 



% — 



2 2A + 1 



i 2A+1 / \ ■ \ 3 

- E (4.0)-i-w + ^l) 

j=0 v 7 
^ 2A+1 

+ 6 E Hmax(j)+Cin(j)-2) (B.14) 



where 



Cin(i) = - o-j+i + 3/2, i" ax (j) = //j- + 0J - 1/2 (B.15) 

are, respectively, defined as the minimum and maximum of i satisfying (i, fa) G 
D{ji) /C*2(/i) in [1 + fjfj]. In the last equality in ( IB. 141) . we have used ( IB. 121) . 
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Substituting (lB~T5l) into (Elij) and using (191 p . S'^_-^k) is expressed as 

. 2A+1 

= a E [(# + *i " V(2A + I)) 3 - K - ^') 3 " 2 ^ 



+ 



M(M 2 - 1; 



A 2 M 1 
+ 



3 (2A + 1) 2 6 

Substituting (IB. 131) and (IB. 161) into (1B.6j) . we obtain 



2A+1 



= -(2A + 1) 2^ I T — 

i=i V 



3L 2 



+ 



2tt 2 (2A + 1) 2 M(M 2 - 1) 2vr 2 A 2 M 



3L 2 



L 2 



— V 

=£jv(g) 



from which and (|63|) . the expressions (1921) and (1931) follow. 



(B.16) 



(B.17) 



Appendix B.2. Matrix elements 

Here we describe X^, y M (r), Z M (r) in the matrix elements in terms of the renormalized 
momenta (|9~T|) and the spin variables (|89|) and (1901) . First we consider X M . The product 
with respect to s — (i,j) G -D(/i) \ C^A*) is taken within each column, and then taken 
over the column. In the j-th column, this condition is equivalent to even (odd) i when 
j is odd (even). The square s = (fi'j,j) belongs to D{ji) \ (^(/i) when j 6 P and 
s = (fij — l,j) belongs to D(/j) \ C^/i) when j e Q. The maximum value of i in the 
j-th column is thus expressed as 

'max = l^j ~ ^o-j<r 2 \ + 2- (B.18) 

The contribution to X M from j-th column is then given by 

'max imax/2 



U (-«(i - 1) + i/2) = (i'-aO'-l)) 



i'=l 



for odd j and 



n (-«(j - 1) + 1/2) 



i:odd 



r[w/2 + i-q(j- i)] 

T[l-a(j-l)] 

(imax+l)/2 

(i' - 1/2 - a(j - 1)) 

i'=l 

r[w/2 + l-a(j-l)] 



(B.19) 



m/2 --a-i)] (R20) 

for even j. We introduce a dummy index i' = i/2 in ( IB. 191) and z' = (z + l)/2 in ( IB. 201) . 
respectively. From (1B.19j) and (1B.20j) and with use of 

r[l-a(j-l)] n r[l/2 - a(j - 1)] 

je[l,3,-,2A+l] je[2,4,-,2A] 
2A+1 

= nr[j7(2A + l)], (B.21) 
i=i 
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we obtain 

= l\ ry/(2A + 1)] ' (a22) 

With use of flB~T8]) and (lB~22l becomes (|9l. 

Similarly, we can write 3^(r) in terms of jjfj and Oj for j G [0, 2A + 1], as shown 
below. When G C 2 (fi), both i and j are odd or even. The maximum value i max 
of i in j-th column is given by 2 max = fi'j — S ajCT0 . The contribution to Y^{r) from j-th 
column is given by 

imax («max + l)/2 

I I Hj - 1) + r + (N - i)/2) = H {a{j - 1) + r + (N + l)/2 - i') 



i:odd i'=l 

r[c*0'-l)+r + (iV+l)/2] 



T[a(j -l) + r + (N- z max )/2] 



(B.23) 



for odd j and 



imax imax/2 

Y[( a {j-l)+r + (N-i)/2)= J] (a(j - 1) + r + JV/2 - i') 

i:even i'=l 

T[a(j -l)+r + N/2] 



T[a(j-l)+r + (N- w)/2] 



(B.24) 



for even j. With use of (191]) and the relation 

J] r[a(j-l)+r+(N + l)/2) J] F[a(j - 1) +r + N/2) 

je[l,3,-»2A+l] je[2,4,-,2A] 
2A+1 

= Y[T[a + r + N/2 + (j - 1)/(2A + 1)], (B.25) 
j=i 

the expression Y^{r) is obtained as ([93]) . 

Next we consider the expression for Z^ir) in terms of {jl'j,aj}. It is convenient to 
decompose -D(/i) as shown in Figure |BT1 Correspondingly, Z^ir) is rewritten as 

2A+1 2A+1 

z,(r) = n n n + r + 

j=l k=j s£H 2 {^)nD jk 

with 

D jk ={s = G [1 + /4 +1 ,/4]} . (B.27) 

Within Djk, the arm length a(s) is constant (= k — j) and hence the squares which 
belong to H 2 (fi) and -D(/x) \ H 2 (fi) are, respectively, aligned alternately. Let 2 m i n and 
Vax be, respectively, the minimum and maximum values of % in s G H 2 (fi) D -Djfc. 
Thus the product over s = (i, j) G H 2 (fi) fl .D^ is expressed as the product over 
i' — 0, 1, • • • , (i max — Vin)/2 with i = 2%' + 2 min . The contribution to ^(r) from Dj k is 
given by 

i [ (aa(s)+r + l(s)/2) 

sGH 2 (fJ,)nD jk 
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B 



Figure Bl. Young diagram of a quasi-hole state for A = 1 is decomposed into Dj k 
for calculation of 



t^max ^min)/2 

= ^ (a(k-j)+r + /4/2 - i min /2 - «') 

i'=0 

= r[q(fc - j) + r + /^/2 - i min /2 + 1] 
r[a(fc-j)+r- + ^/2-w/2] 
The maximum value z max of i in H 2 (fi) fl Pjfc is expressed as 



(B.28) 



«max = /4 - (B.29) 

as shown below. 

For s G Djfe, a(s) is — j and Z(s) = fi'j — i and hence 7t(s) = a(s) + l(s) + 1 is 
written as 

h(s) — f/j— j + k — — fjfj— j — (/jf k — k) + fx' k — i + 1. (B.30) 

When (j, k) G (P, P) or (Q, Q), both //• — j and n' k — k are odd or even and hence for 
G H 2 (fi)nD jk , 

H' k -l = i mod 2. (B.31) 

When (j, k) G (P, Q) or (Q, P), on the other hand, the sum of ^ — j and \j! k — k is odd 
and hence, 

fif k = i mod 2. (B.32) 

The two relations (IB. 311) and ( IB. 321) are expressed as ( 1B.29I) in a unified way. In a 
similar way, i min is expressed as 

Vin = /4+i + 1 + <W fc+ i ( B - 33 ) 

respectively. 
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The expression f]B.28|) is further written as 
\\ (aa(s)+r + l(s)/2) 



s&H 2 ([t)nD jk 

r[(g - K+x - <WJ/ 2 + r - <* + V2] 



(B.34) 



(CI) 



r[(#-ft + ^)/2+r] 
with use of (!R29l) . fjR33l and (JSU). From (jR34p , the expression ([961) for Z^r) follows. 

Appendix C. Spectral Weight 

The triple integral in the spectral function for A = 1, 

dui J du 2 J du 3 5 I p - — (ui + u 2 + u 3 ] 

x s(e-^^- (3 - uj - u 2 2 - uj)^J F(u u u 2 ,u 3 ) 
with 

F(u 1 ,u 2 ,u 3 ) 

= \u 2 - n 3 | 4/3 It*! - u 3 \~^ \ Ul - u 2 \~ 2 " f[ (1 - u>y l/3 (C.2) 

i=i 

reduces to an integral on a curve determined by the sphere, plane, and the cube 

e = ^(3-ut-ul-ul) (C.3) 

P=^-(u 1 + u 2 + u 3 ) (C.4) 
K|<1, i = 1,2,3. (C.5) 



The cross section between ( 1C4[) and (1C 5|) is triangular when 



nd/2 <p< 3nd/2 (C.6) 

and hexagon when 

< p < nd/2, (C.7) 

as shown in Figure ICU 

In the following, we consider the case (lC.6j) only. The case (IC.7I) can be discussed 
similarly. The integral in (jC.lft is on a circle when 

7cd 3nd 3 ( 7rd\ 2 3, ,, 2 ^ ^ 2 9/ ,\2 n \ 

Y < P < and - - U - —J + -(vrrf) 2 < e < -p 2 + -(nd) 2 , (C.8) 

and it is on three disconnected pieces of arc when 

ird/2 <p< 3ixd/2 

and 

"3(P - -rf) 2 + ^) 2 < e < -\ (p - ?f) 2 + 3 -(nd) 2 , 
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Figure CI. Cross section of the 
plane (|C.4|) and the cube (|C.5|) . 
The cross section is triangular 
PQR when p S [nd/2, 3i:d/2] and 
hexagon when p £ [0, nd/2]. 



Figure C2. Cross section of the 
plane (|C.4|) and the cube (|C.5jl for 
p e [nd/2,3ird/2}. 



as shown by bold curves in Figure IC2l The integrand F(ui, w 2 , u 3 ) in (1C.1I) diverges 
when 

ui = u 3 , Ui = u 2 or \uj\ = 1 for j = 1, 2, 3, (C.9) 

which are represented by solid lines in Figure IC2I When the contour of the integral 

passes near the crossing points of the above lines ( 1C.9I) . the spectral function becomes 
singular. Those crossing points are given by 

j 2p 2p 2p \ 

1,1.-2 + 1), (C.ll) 

1P 1 +^,ll (C.12) 



2 vrrf' 2 vrrf 

and their equivalent points. They are depicted in Figure IC2l by dots. The singularity of 
the spectral function near the upper edge e = (37r<i/2) 2 (i — (2p/(ird)) 2 ) of the support 
comes from fIC.lOj) . The point fIC.llj) yields the singularity of A(e,p) near the lower 
edge e = 3(7rd/2) 2 (2p/(7rd) — l)(3 — 2p/(nd))). flC.12p is relevant to the singularity near 



2 V 2 J \ 2 

The singularities near the upper edge and lower edge have been obtained in earlier 
papers[B21 US]- We thus consider the singularity when Se = e — e p ~ in the following. 
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Figure C3. Cartesian coordinate {e^e^e^} and cylindrical coordinate {e p ,eg,e^}. 

Now we introduce the Cartesian coordinate ■ ej = Sij i,j = 1, 2, 3 and define the 
vector 

u = u x ei + u 2 e 2 + u 3 e 3 . (C.14) 
We define another Cartesian coordinate 

eg = -e 1 /V6 - e 2 /V~6 + 2e 3 / \/6, 
e v = -ei/y/2 + e 2 / V2, 

e c =(e 1 + e 2 + e 3 )/V3 (C.15) 
and circular coordinate on eg — e v plane 
e. p = eg cos 9 + e v sin 6, 

e e = — eg sin 6 + e v cos 6. (C.16) 

Note that {eg, e n , eg} forms a cylindrical coordinate as shown in Figure IC31 In terms 
of (IC.15I) and flC.16j) . we rewrite u as 

u = u p e p + wgeg. (C.17) 



In the following, we change the variables of integral in flC.lj) from (ui,u 2 ,u 3 ) to 



(u p ,9,u c ). From (!CTl4l) -( !CTT7l) . we obtain 



Ui — u ■ e\ — u p (— cos 8/V6 — sin 8/ V2) + u^/V3 
u 2 = u ■ e 2 = u p (— cos 6/VQ + sin 0/y/2) + u c /y/3 

u 3 = u- e 3 = 2u p cos e/Ve + ujVz. (c.i8) 

The spectral function ( 1C.1I) is rewritten as 

\^3ndut 



A (e,p) = C J dupUp J d9 J d«g 5 — 



x 5(e-3(ird/2) 2 (3-uJ-uJ)) 

xF(ui,u 2 ,U3)- (C19) 
From the delta functions, u p and «g are forced to be 

u p = v/3-4(e + p 2 )/(37r 2 d 2 ), m c = 2p/( v / 37rrf), (C.20) 
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respectively. As a result, A (e,p) becomes 

A-, N AC 

A (e,p) 



d0F(«i(0),U 2 (0),U 3 (0)) 



(C.21) 



with 



(0) = v/(l/2-j9'/3) 2 -2(5eV9(- cos 9 - V3 sin 0) + 2p'/3 



u 2 (9) = Vi 1 / 2 - P'/ 3 ) 2 ~ 25e'/9(- cosfl + v 7 ^ sin 0) + 2p'/3 
«3(0) = - 2p'/3) 2 - A5e'/9 cos + 2p/3. (C.22) 
Here we have introduced Se' = Se/ijid) 2 and p' = p/(7rd). The integral in ( 1C.21I) runs 



over 



0G[O,2tt], for 5e>0 

G [0 O , 2tt/3 - 0„] U [2tt/3 + 0„, 4vr/3 - 0„] U [4vr/3 + 0„, 2vr - O ] (C.23) 
for <5e < with 



9 = arcsin 



The point (IC.12I) and equivalent points 

1 p 1 p 

2 Vd"~2 %d 



-2Se 



(3/2 -p') 



25e 



1 



1 p 

2 ^ 



, 1, 



1 p 

2 Vd 



(C.24) 



correspond to <5e' = and 9 = 0, 27r/3, 47r/3, respectively. The most singular 
contribution to A~(e,p) for 5e = comes from the vicinity of 9 = 0, where the factors 
| Hi — W2|~ 2 / 3 (l — M3) 173 become singular in F(ui,U2,u^). Therefore, we approximate 
F(u l} u 2 ,u 3 ) as 

Isin0|- 2 / 3 



Foc(w 1 (0)-« 2 (0))- 2 / 3 (l-w 3 (0))- 1/3 



ex 



[1 -gcos9)V 3 



1 — Se'/(p' — 3/2) 2 . A (e,p) is evaluated as 
Isin0|- 2 / 3 



with g = 

A~(e,p) oc 
for Se > and 5e ~ and 
A'(e,p) oc 



[l-gcos9) 1 / 3 

1 sm^r 2 / 3 



d0 



d0 



(C.25) 



(C.26) 



(C.27) 



(1 -#cos0) 1 /3 

for 5e < and 5e ~ 0. Here C is a cut-off angle of the order unity. We can take 
9 C = 7r/3, or 7r/2, for example. First we evaluate ( 1C.26I) . Introducing t = tan(0/2), 
t c = tan(0 c /2) and g = (1 - #)/(l + (10261) becomes 

to r 2/3 (1 + t 2N-4/3 



A (e,p) oc 



(* 2 + s) 1/3 



-df. 



We rewrite fl028l) as 

£ -2/3 (l + t 2N-4/3 



(t 2 +^)V3 

(* 2 + s) 1/3 



df 



df + 



fc r 2 / 3 ((i + t 2 )- 4 / 3 -i) 



df 



(C.28) 

(C.29) 
(C.30) 
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The second integral of the right-hand side converges when g = 0. The first term of the 
right-hand side, on the other hand, is rewritten as 

"t c g~ 1/2 T -2/3 poo --2/3 

«- 1/6 / VTW dT ~ rl ' e l (r 5 TT) I 73 dT ' (a31) 

Since g oc Se, we arrive at 

A-{e,p) oc{5e)- 1/6 (C.32) 

when 5e > and 5e ~ 0. 

The singularity near <5e = and 5e < is evaluated similarly. Introducing 
t = tan(0 o / 2 ), (E22l) becomes 

When g ~ 0, equivalently 5e ~ 0, tog' 1 ' 2 is the order of unity and hence, 

A-(e,p) oc (-5e)- 1/6 (C.34) 
when 5e < and 5e ~ 0. 
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